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Optimum Baffle 


Baffle (Ht/Wth) = 0.83 
Percent of Combustor 


Diameter 


=35.5 




e3 a> 

> C 

g '2 

Ut nj 

H # g 

^ *53d 

fl 

2 W 

cs __ 

co 13 

'■o .a 


& 


a> 


-g W 

$5 cm 

® f 

j: -g 

O § 

e a 

« t; 
£ as 
x o, 
« i> 
CQ Q 


•a ^ 

& 3 

> Tf 

'£ ^ 
P .2 

0> 5 

3 ® 

CO H3 

T3 5 
a JS 

cs a> 
*3 £ 

iu C 

0 


8th International Conference on Mathematical and 

Computational Modeling 






3 

O V 
r> ca 

3 g 

« ’f 

2 S 

3 

O 

43 

CA 

<X> 

N 

*CA 

s 

«t-n 

<3 
42 

a „ 

3 -a 

CA 
<V 

43 

CJ 

fl 


m 

68 

U 

a 

£ 

OD 

• pH 

0 ) 

p 3 


a 

a 

o 

4 > 

43 

H 


CA 

0 > 


a s 

*> ‘a 

fS pQ 


o> 

CA 

3 

S*-H «0 

O Oh 

si 

0 ) = 

§ a 

&£ 

a w 

2 £ 

h 42 

3 ft 

Si 


<a> 

CA 

03 

43 

Oh 

i 

«*h 

0 

1 

3 

_© 

4 ) 

43 

H 


CA 
4> 

'V 

4> 

!- 

O 

a 

CA 

CA 

4> 

§S 


3 

40 


3 

&> 

a 

<v 

00 

3 

3 

U 

U 

3 



FLOW VISUALIZATION AND COMPUTATIONAL 
STUDIES OF A REVERSE FLOW 
COMBUSTOR 


B. Ghorashi, F. Reardon, G. McBeath 

and 

K. Chun 

Chemical Engineering Department 
Cleveland State University 
Cleveland, Ohio 44115 


32nd Heat Transfer and Fluid Mechanics Institute 





nonetheless, additional optimization 
studies are needed 



proceedings of the 

32nd HEAT TRANSFER AND 
FLUID MECHANICS INSTITUTE 

Held at 

California State University, Sacramento 
June 6 & 7, 1991 


Edited by 

Frederick H. Reardon 
and 

Ngo Dinh Thinh 

Department of Mechanical Engineering 


School of Engineering and Computer Science 
California State University, Sacramento 
Sacramento, California 



CONCLUSIONS 


PATTERN DEPICTED BY SMOKE VISUALIZATION 
WAS DESIRABLE FOR RAPID MIXING OF FUEL AND 
AIR 


MORE EFFECTIVE ATOMIZATION OF FUEL IS 
REQUIRED 


LOW INJECTION VELOCITY AND LARGE SPRAY 
ANGLE PRODUCED A MORE UNIFORM 
TEMPERATURE DISTRIBUTION AT EXIT 



FLOW VISUALIZATION AND COMPUTATIONAL STUDIES 
OF A REVERSE FLOW CIRCULAR COMBUSTOR 

B. Chornhi, F. Reardon*, C. McBeath. and K. Chun** 

Chemical Engineering Department 
Cleveland State University 
Cleveland* Ohio 44115 

ABSTRACT 

Computational Fluid Mechanics Studies using the Fluent computer program and 
flow vis ualiza tion studies were undertaken to support the development of a reverse- 
flow circular combustor design. The reverse-flow drcular combustor combing a 
drcular geometry with multiple fuel injection nozzles in Mini -Combustion Zones (MGZ). 
An overall combustor model was developed and used to determine the effects of. fuel 
injector design on the temperature distribution throughout the combustor, with 
particular attention to the combustor exit Low injection velodty and Urge injection 
spray angle were found to produce a more uniform temperature distribution at the exit 


The flow pattern depicted by the smoke visualization technique showed a 
behavior which is very desirable for the rapid mixing of fuel and air. The atomization 
of fuel was visualized by injection of a yellow dye smoke in the liquid side of the air- 
blast injector. The portion of the combustor volume that was investigated consisted of 
four baffles on the outermost drcular surface and two baffles on the innermost surface. 
The flow exhibited a rapid circulation pattern inside the MCZs. The channel flow also 
exhibited a high degree of interaction with the MCZs. Experimental tovestigation of the 
flow pattern in the combustor showed similarities to the computational fluid dynamic 
results. 


INTRODUCTION 

The development of efficient air-breathing jet propulsion engi nes t hat produce 
low levels of pollutants in the exhaust will require an improved understanding of the 
flow and reaction processes that take place in combustion chambers. These processes 
indude atomization and vaporization, when Kquid fuels are used, and transport mixing, 
and chemical reaction, for both liquid and gaseous fuels. Computational fluid mechanics 
£s a technique for studying these processes without the need for excessive experimental 
runs. 


A flow visualization study of the mixing pattern in a reverse flow combustor 
was conducted (Figure 1). A transparent model was constructed to examine the flow 
pattern and mixing of side- mounted turbulent jets and their interaction with the main 
channel flow. The portion of the combustor volume that was investigated consisted of 

• Mechanical Engineering Department California State University, Sacramento 

- NASA, Lewis Research Center 
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SCALAR COMPUTATION MODULE FOR CFD CODES WITH FLUID AND KINETIC EFFECTS 
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Cleveland, Ohio 44115 

and 

Kue Chun 
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ABSTRACT 


A model was developed to interface with a CFD, k-c based code. 

A converged solution from the CFD code is the input to the post- 
processing model for prediction of thermal N0 X . The model uses a # 
decoupled analysis to estimate the equilibrium level of (N0 x ) e which 
in turn is used to predict its rate of formation, (dNO x /dt) e , at each 
node. The rate is fixed in the N0 X production rate plot by estimating 
the time to reach equilibrium, (t) e q, by a differential analysis based 
on the reaction: 0 + N 2 NO + N. The rate is integrated in the non- 
equilibrium time space based on the residence time at each node in the 
computational domain. The sum of all nodal predictions yields the 
total N0 X level from the combustor. Advantages of this model include: 
(1) use of all parametric variables of interest for optimization (T,P, 
N 2l 02,N0 ef /i,p,i:i ce ,V m ,d n ,N ref D e f ,A 0 ,A) , (2) turbulent mixing effects, 
(3) geometric effects, (4) lean-to-rich burn combustion conditions, 
and (5) the potential to work with most k-c based CFD codes. 



A COMPARISON STUDY OF STRAIGHT AND CIRCULAR TUBS 

COMBUSTORS 


Bahman Ghorashi, Sastry Taruvai and Georgio McBeath 
Department of Chemical Engineering 
Cleveland State University 
Cleveland, Ohio 44115 


and 

Kue Chun 

NASA Lewis Research Center 
Cleveland, Ohio 44102 


ABSTRACT 


A Numerical study was conducted using a CFD, k-e based code [ ] . 
The objective of the study was to compare two combustor 
configurations i.e., a straight tube combustor and a reverse-low 
circular combustor. Each of these geometries were further 
subdivided into one containing the concept of Mini-Combustion- Zones 
(MCZ) and the other without any compartmentization. The model with 
MCZ's contained internal baffles and multiple fuel injection ports. 
The other model contained only a series of multiple injection 
ports. A one-to-one comparison was made between straight and 
circular tube combustors with and without MCZ's. The objective 
was to identify a geometry that could provide complete and 
immediate combustion of fuel by using the entire combustor space as 
well as a uniform temperature at the exit of the combustor. The 
uniformity of temperature was used as a measure of the goodness of 
mixing . 


INTRODUCTION 



A COMPARISON OF FLOW VISUALIZATION EXPERIMENTS WITH 
COMPUTATIONAL FLUID DYNAMIC SIMULATIONS OF A CIRCUIAR 
COMBUSTOR WITH MINI -COMBUSTION ZONES 

Giorgio McBeath, Bahman Ghorashi and S as try Taruvai 
Cleveland State University 
Cleveland, Ohio 44115 U.S.A. 

and 

Kue Chun 

NASA Lewis Research Center 
Cleveland, Ohio 44135 U.S.A. 

ABSTRACT 

Flow visualization experiments were performed to study the nature 
of complex flow pattern in a circular combustor containing mini- 
combustion zones (MCZs) . The MCZs were structured by bell-shaped 
baffles located along the inner walls. The premise was that a series 
of baffles placed along the wall of the combustor could enhance mixing 
and better utilize the combustor volume. Complex flow patterns were 
visualized and were photographed with a high speed camera. The 
results were analyzed and compared to computational fluid dynamic 
predictions from the Fluent code [1]. The simulations were performed 
at idle and cruise conditions. 

The flow visualization studies examined the mixing effects 
of two air-blast nozzle designs. One nozzle was a standard expanding 
turbulent jet atomizer and the other was an air-blast atomizer for 
turbine engines with a 45 degree swirl vane generating 20 to 60 micron 
droplets. Using the uniformity of temperature as an indirect measure 
of goodness of mixing, it was determined that a combustor with mini- 
combustion zones has definite advantages over a straight-through 
combustor. Idle and cruise conditions showed similar results in flow 
pattern in the range of M - 0.1 to 1.0 as determined by simulations. 

This study also examines the interaction of the airblast nozzle 
in the mini -combust ion zone (MCZ) with the main channel flow. Rapid 
mixing of combustion products with the main channel flow and the 
minimization of temperature non- uniformity would minimize NOx 
generation. The results show a mechanism of mixing which is conducive 
to reducing nitric oxides in turbine combustors. 



OPTIMIZATION STUDIES OF A REVERSE FLOW CIRCULAR 
COMBUSTOR WITH MCZ'S 


Bahman Ghorashi, S as try Taruvai and Georgio McBeath 
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Cleveland, Ohio 44115 


and 


Kue Chun 
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ABSTRACT 


A numerical study was conducted using a CFD, k-e based code [ ] . 
The objective of the study was to optimize the design of a reverse 
flow circular combustor [ ] . The circular combustor combines a 
circular geometry with multiple fuel injection nozzles in the Mini- 
Combustion-Zones (MCZ) . The concept of MCZ's was incorporated by 
having a series of injection ports and baffles, figure [ ] , whereby 
each injection nozzle and its surrounding baffles would comprise an 
MCZ. 

An overall combustor model was developed and was used to determine 
the effects of nozzle position, fuel spray angle, baffle positions 
with respect to each other, number of baffles, dimensions of 
baffles and the location of baffles. Initially the baffles were 
placed on top and bottom surfaces of the combustor and the 
optimization studies were performed ( ] . Industrial constraints 
were later introduced into the design and further studies were 
performed with baffles placed only on one side. The objective of 
this study was then to compare these two final design 
configurations and determine which of these geometries could 
provide the requirements such as low pressure drop within the 
combustor, uniform temperature at the exit of the combustor, (a 
measure of goodness of mixing) , as well as complete and immediate 
combustion of fuel by using the entire combustor space. 



A NUMERICAL SIMULATION OF THE EFFECTS OF POSITION AND 
SPRAY ANGLE OF THE FUEL INJECTION NOZZLE 
ON MIXING IN A CIRCULAR COMBUSTOR 

Bahman Ghorashi and S as try Taruvai 
Chemical Engineering Department 
Cleveland State University 
Cleveland, Ohio 44115 

ABSTRACT 

A Numerical study was conducted using Fluent computer program developed by Creare, Inc. 
The objective was to determine the flow pattern between two sets of baffles in an experimental 
reverse-flow circular combustor, designed for enhanced mixing and NO x reduction,by varying the 
position and spray angle of the fuel injection nozzles. The segment of the combustor that was 
investigated consisted of two baffles positioned on the upperface and two on the Iowerface of the 
circular combustor. Heated air at 300°F and 3 atmospheres pressure was introduced into the 
combustor with an inlet velocity of 100 ft/sec. The fuel, pentane, was introduced into the main 
stream via fuel injection nozzles. There were four nozzles, two on the upper and two on the 



NUMERICAL SIMULATION OP FLOW AND THE EFFECT OF BAFFLE 
ARRANGEMENT ON PRESSURE DROP AND TEMPERATURE PATTERN 

IN A CIRCULAR COMBUSTOR 


Bahinan Ghorashi and Sastry Taruvai 


Department of Chemical Engineering 
Cleveland State University 
Cleveland, Ohio 44115 


ABSTRACT 

A numerical simulation study was conducted to determine the effect 
of baffle arrangement on pressure drop in an experimental reverse- 
flow circular combustor with several Mini-combustion Zones (MCZ) . 
The Fluent computer code developed by Creare Inc. , was' used for 
this study. Three parameters which appeared to significantly 
affect the pressure drop were varied, and pressure drop for each 
case was numerically determined. The portion of the combustor 
volume that was investigated had two baffles on the upper face and 
two baffles on the lower face. Heated air at 1000°F with an inlet 
air velocity of 100 ft/sec and 1 atmosphere pressure was introduced 
into the combustor. Fuel (pentane) at 80°F was introduced via 
nozzles which were placed between the baffles. The following 
changes were then made in the baffle arrangement and the effects on 
the pressure drop and temperature distribution were evaluated for 
each of the following cases: 

1. The positions of the baffles on the upper and lower face were 
changed . 

2. The dimensions of the baffles were varied. 

3. The shapes of the baffles were altered. 

The effects of the above mentioned design alterations are 
summarized. 

KEYWORDS 

Mini-combustion-zones; baffles; pressure drop; combustor. 
INTRODUCTION 

A reverse-flow circular combustor based on the design proposed by 
Ghorashi [1,2,3) in 1987 was studied. This approach introduces the 
concept of Mini-Combustion-Zones (MCZ) to use the entire reactor 
volume via radial injection of fuel along the axial main flow. In 
addition to this, the use of Mini-Combustion-Zones generates 
recirculation within the MCZs, making each MCZ to behave like a 
well-stirred reactor. *- 
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Summary of LDV Results 

Main flow and nozzle flow mixing is indicated 
by high turbulence intensity values (peaks) at 
the same location. 


Higher velocities were evident near the inner 
radius of the combustor. 


Multiple high turbulent intensity peaks in the 
MCZ indicate significant mixing is occurring. 


Shear layers from baffle peaks interact with 
nozzle flow increasing local mixing 


Mainflow air adds to the uniformity of mixing 
across the Y-cross section and enhances the 
overall mixing in the Y-plane (Area under 
curve). 


Main flow pressure-increases result in higher 
axial velocity and lower radial velocities. 
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COMPARJSION OF COMPUTER SIMULATION RESULTS 
FOR C«H 1 8 MASS FRACTION 
MASS FRACTION DIMENSIONLESS 
ORIENTATION X, PALNE 6, 3 - D DOMAIN 

(a) - Circular, no Baffles (b) - Circular, Baffles 
(ci - Straight, no Baffles (d) - Straight, Baffles 
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COMPARISON OF COMPUTER SIMULATION RESULTS 
FOR TEMPERATURE PATTERN 
TEMPERATURE IN DEGREE RANKINE 
ORIENTATION X. PALNE 6, 3 - D DOMAIN 


(a) - Circular, no Baffles (b) - Circular, Baffles 
(cl - Straight, no Baffles (d) - Straight. Baffles 
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COMPARISON OF COMPUTER SIMULATION RESULTS 
FOR C02 MASS FRACTION 
MASS FRACTION DIMENSIONLESS 
ORIENTATION X, PALNE 6. 3 - D DOMAIN 

(a) - Circular, no Baffles (b) - Circular. Baffles 
(c) - Straight, no Baffles (d) - Straight. Baffles 
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THERMAL N0 X PREDICTION MODEL 


PROCEDURE: 

Converged Fluent Case (k-epsilon Model) 

1. INPUT: 

Pressure 
Temperature 
Oxygen Concentration 
Nitrogen Concentration 
Eddy Dissipation 
Effective Viscocity 
Density 

Velocity Magnitude 



2. Calculate Effective Diameter (D e ) and the Cross- 
Sectional Area (A x ) of the Combustor at Each Node 

Plan 

D e = 4 * Area / Perimeter 
A x = L * W 

3. Calculate Kolmogorov Length Scale at Each Node 

Ao = [(|X/p)3/e]0.25 

4. Calculate Turbulent Reylolds Number at Each Node 

Nfe = [D e |V m |p/(i] 

5. Calculate A.pm at Each Node 

^pm = V(|i*A / D e |V m |p) 

3. Calculate the Residence Time for Each Node 


tRT = d / | V m | 


7. Test for the Degree of Mixing at Each Node 
( if 0 < A.pm/ Xo < 1, go to step A 
if ^pm/^o = 0, go to step B) 


A. Good Mixing = 2.78 * [Hung's Eq.] (Smith, et.al.) 

B. Diffusion Limited Case = 1.00 * [Hung's Eq.] 





8. Correct for Nitrogen Conversion at Each Plane 
(Mass Balance at Each Plane to Determine the 
Nitrogen Consumption) 

Note: Correction for Nitrogen Conversion at each 

Node Might Yield More Accurate Results, 
However, it would be Computationally 
Inefficient. Nontheless, the latter should 
be tested. 


[N 2 J 0 = 0.79 - 1/2 * Z[NO] t - I[N] t = [N 2 ] t 
[NO] t = 2yt - y/a(l - e-at) 

[N] t = 2pyt/a + y/a(l - e-«t) 


where a , p and y are Determined from the 
Arhenious Equation 


9. Calculate the Equilibrium Values of Oxygen, 

Nitrogen and Nitric Oxide with 18-Step Octane 
Combustion Model (Kopa, et.al.) 


hrtma i, «sul id lie aimbm of a'cmonaoeta 

** munbcr of ramie meets. SLa ti. Car elements in 
atomic term: CH.O.K. 
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10. Calculate the Rate of Nitric Oxide Formation at the 
Intersection of Equilibrium NO and NO Rate-Line, 
Using the Hung Equation at Each Node 

cIXno / dt = P 1 / 2 [(aX n 2 X 0 2 1/2 * bX N 0 2 /X o 2 i/2) / (i + cX NO /X 02 )] 
Where a,b,c is determined from the Arhenious Equation 



11. Estimate the Relaxation Time Via the Differential 
Analysis of the Controlling Mechanism 


tcalc - [NO] e /2y 


12. Calculate the Actual Time Based on Two-Step 
Mechanism 

O + N 2 = NO + M 
N + 0 2 = NO + O 
*RX = t caJc /2 


13. Construct Semi-Log NO x Formation Rate Curves 
(refer to Step 10) 
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14. Determine Non-equilibrium NO x Formation Rates 

Note: Interpret Curve to Find Non-equilibrium NO 

Rate at Lower Limit [tRj * 0.01] and at 
Upper Limit [tRj]- 


15. Integrate the Rate with respect to Time 

NO(i,j)= JdNO/dt*dt=0.5’(t R T-0.01*t RT )*AdNO/dt 

16. Sum all the Local NO x Values to Estimate the Total 
Thermal NO x 

m 

NO x = Z NO(i,j) for K = 1 to N 
i.i - 1 




ALTERNATE METHOD 


Procedure 

Converged FLUENT Case 

1 . INPUT 

Pressure 
Temperature 
Oxygen Concentration 
Nitrogen Concentration 
Velocity Magnitude 

Assume Constant Nitrogen Concentration 

2 Use Global Forward-Rate Reaction to 
Estimate the Local NO x Formation 

[NO] = Ke(-k2/T)[N 2 ][0 2 ] 1/2 t 
Where [ ] = mole fraction 

T = Temperature 
t = time 
K,K2 = constants 

3. Integrate with respect to Residence Time 



4. Sum all the Values in order to Obtain the 
Total Thermal NO x 

Limitations of this Method 


1 . Temperature Range 


?«mo«r»4urt. *f 

3092 1272 34S2 3412 3812 



2030 oc < T < 1700 °C 

2 . Assumptions 

dX N0 /dt = Pi/2[(aX N2 X 0 2 1/2 -bX N o 2 /Xo2 1/2 )/(1+ cX N0 /X 0 2 )] 

if Xnq small, then (Xnq) 2 = 0 
if Xno is small, then Xnq/X 0 2 ~ 0 

[NO] = Ke(-k2/T)[N 2 ][Q 2 ]i/2t 


ORIGINAL Paqb is 












The equivalent molecular diameter is given by equation 



The fnass diffusivity can be calculated from: 




The Schmidt number of air is given by 


(Sc) 


air pD 


AB 


(Sc) 


air D 


AB 


For the high-S chmi dt-number cas e, the time constant can_be_ obtained by 




no mix ing in the limit of infinite Schmidt number is now satisfied (r = oo, 
•which gives /, = 1 for all time). 
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The local mixing pr oble m is simplified 


h = 


1 

1 + t'/t 


Isotropic theory as previously presented could be applied to mixing-vessel 
problems with some degree of approximation. However, in this case the turbulent 
field will decay in certain regions of the mixer. Probably the only area which 
would approximate a constant turbulent field is that in the immediate vicinity 
of the stirrer. 
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